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We find that the statistics of levels undergoing metal-insulator transition in systems with Gaussian 
disorder and non-interacting electrons behaves in a way similar to that of the single parametric 
Brownian ensembles [1]. The latter appear during a Poisson — * Wigner-Dyson transition, driven 
by a random perturbation. The analogy provides the analytical evidence for the single parameter 
^ scaling of the level-correlations in disordered systems as well as a tool to obtain them at the critical 

• point for a wide range of disorders. 
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The spectral correlations of a disordered system are very sensitive to the behavior of its eigenfunctions. The 
presence of disorder may cause localized waves in the system, implying lack of interaction between certain parts. This 
is reflected in the structure of the Hamiltonian matrix which is sparse in the site representation. The degree of sparsity 
of the matrix is governed by various system parameters e.g. dimensionality, shape, size and boundary conditions of 
the system. The variation of the disorder-strength can lead to a metal- insulator transition (MIT), with eigenfunctions 
changing from a fully extended state (metal) to a strongly localized one (insulator) with partial localization in the 
critical region. The associated Hamiltonian also undergoes a transition, (in effect only due to variation of the relative 
i ^ i ■ strength of its elements) , from a full matrix to a sparse or banded form and finally to a diagonal matrix. The statistical 
studies of levels for various degrees and types of disorders as well as system conditions require, therefore, analysis of 
different ensembles. Here the nature of the localization and its strength is reflected in the measure and the sparsity 
of the ensemble, respectively. Our objective in this paper is to obtain a mathematical formulation for the level- 
■ correlations, common to a large class of system conditions (with Gaussian type randomness); the system information 
7-H ', enters in the formulation through a parameter, basically a function of various system parameters influencing the 
t-H ■ localization. 

Recently it was shown that the eigenvalue distributions of various ensembles, with a multi-parametric Gaussian 
measure and independent matrix elements, appear as non-equilibrium stages of a Brownian type diffusion process [2]. 
Here the eigenvalues evolve with respect to a single parameter which is a function of the distribution parameters of 
the ensemble. The parameter is therefore related to the complexity of the system represented by the ensemble and can 
be termed as the "complexity" parameter. The solution of the diffusion equation for a given value of the complexity 
parameter gives the distribution of the eigenvalues, and thereby their correlations, for the corresponding system. A 
similar diffusion equation is known to govern the evolution of the eigenvalues of Brownian Ensembles (BE) [1,3] and 
many of its solutions for various initial conditions have already been obtained [4] . The analogy can then be used to 
obtain the level-correlations for the Gaussian random matrix models of the disordered systems with non-interacting 
electrons. The presence of interactions introduces a correlation between matrix elements of the ensemble representing 
the system; the details of this case are discussed elsewhere [5]. 

The correlations in the single electron spectra of disordered metals are governed by a variety of parameters e.g the 
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associated energy-ranges, degree of disorder, the dimensionality of the system etc. Here the two energy-scales, playing 
the dominant role, are the Thouless energy E c and the mean level spacing A. The E c is given by the time-scale 
needed by the wave-packet to diffuse through the sample. In the diffusive (metallic) regime and for energy-scales 
8E smaller then E c , the spectral correlations are well-modeled by Wigner-Dyson (WD) ensembles [3]. For 5E > E c , 
the statistics deviate from the Wigner-Dyson case, however the deviations are negligible for sample size L — > oo. 
In the localized (insulator) phase too, the correlations are energy-dependent but, in the limit L — > oo, the levels 
are completely uncorrelated and their statistics can be modeled by the Poisson ensemble. However the statistics in 
the critical region near the metal-insulator transition (Anderson type) is different from both Wigner-Dyson as well as 
Poisson statistics and depends on various system dependent features [7,8]. Our study shows that the multi-parametric 
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level-statistics in the critical region can be well-modeled by the single-parametric Brownian Ensembles. 

The paper is organized as follows. The section I contains a brief description of the simplest model of a disordered 
system using independent electron approximation and the equation governing the evolution of its eigenvalues due to 
change of disorder etc. The properties of the BEs useful for present study are given in section II. The section III 
deals with the determination of the single parameter A governing the level-statistics during MIT using BE-analogy. 
It also provides an explanation, in terms of A, of some of the observed features of the AE-statistics. In section IV, the 
AE-BE analogy is used to to obtain the analytical formulation of some of the unknown spectral fluctuations during 
MIT. The section V contains the details of the numerical comparison of the level-statistics of Anderson Hamiltonian 
with that of BEs and reconfirms our analytical results. The studies, during last decade, indicate the surprising success 
of power law random banded matrices (ID system) as a model for Anderson ensembles [6]; as discussed in section VI, 
the success can be explained within A formulation of the level-statistics. 

I. THE MULTI-DIMENSIONAL ANDERSON HAMILTONIAN 

The Anderson model for a disordered system is described by a (i-dimensional disordered lattice, of size L, with 
a Hamiltonian H = J2n enCL n a n — J2 n ^m^ mn ( a n a m + a « a m) m tight-binding approximation. The site energies e n , 
measured in units of the overlap integral between adjacent sites, correspond to the random potential. The hopping is 
assumed to connect only the z nearest-neighbors (referred by m) of each site. In the site representation, H turns out 
to be a sparse matrix of size N = L d with diagonal matrix elements as the site energies Hkk = £fe- The off-diagonals 
Hki describe the interaction between two sites k and I; here Hki for two sites connected by hopping will be referred as 
hopping off-diagonal and the rest as non-hopping off-diagonals. The level-statistics of H can therefore be studied by 
analyzing the properties of an ensemble of (i) sparse real symmetric matrices, in presence of a time-reversal symmetry 
and (ii) sparse complex Hermitian matrices in absence of a time-reversal. 

We consider an ensemble of Anderson Hamiltonians (later referred as Anderson ensemble) with a Gaussian type 
disorder. The site-energies Hkk = are thus independent Gaussian distributions Pkk(Hkk) = e _ ( fffcfc ~ bfcfc ' l 2hkk 
with variance hkk and mean bkk- The hopping can be chosen to be isotropic or anisotropic, non-random or random 
(Gaussian). A general form of the probability density p(H) = Yik l-k<l Pkl(Hkl) of the ensemble, including all the 
above possibilities, can therefore be given by 
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p{H, h, b) = Ccxp[- ]T J2( 1 / 2h ^)(Hkl;s - bkus) 2 ] (1) 

s=l k<l 

with subscript "s" of a variable referring to its components, (3 as their total number ((3 = 1 for real variable, (3 = 2 
for the complex one), C as the normalization constant, h as the set of the variances hki- s =< H^is > an( ^ & as the set 
of all mean values < Hki- S >= bki- s - As obvious, in the limit hki-i,hki-2 — * 0, eq.(l) corresponds to the non-random 
nature of Hki (that is, pki(Hki) = S(Hki — bki))- Note although the non- hopping off-diagonals in Anderson matrix 
always remain zero but the effective sparsity of the matrix changes due to change in relative strength of the diagonals 
and the hopping off-diagonals. Thus, in the insulator limit (with almost no overlap between site energies due to 
strong disorder), the matrix behaves effectively as a diagonal one, the diagonals being very large as compared to 
hopping off-diagonals. In the opposite limit of very weak disorder when an average diagonal is nearly of the same 
strength as an average off-diagonal, the statistical behavior of the matrix is same as that of a matrix taken from a 
Wigner-Dyson ensemble [3]. The latter are the basis-invariant Gaussian ensembles of Hermitian type, with a same 
variance for almost all matrix elements. The statistical behavior of levels in the Wigner-Dyson ensembles depends 
only on their symmetry class and is therefore universal in nature. The three main universality classes are described 
by a parameter (3, basically a measure of the degree of level-repulsion [3]: (i) GOE with (3=1, corresponding to 
time-reversal symmetry and integer angular momentum, (ii) GUE with (3 — 2 and no time-reversal symmetry, (iii) 
GSE with (3 = 4 and time-reversal symmetry but half integer angular momentum. 

A variation of disorder and hopping rate changes the distribution parameters of the probability density p(H) and 
thereby its statistical properties. Using Gaussian nature of p, it is easy to verify that under a change of parameters 
hki — * hki + Shki and bki — > bki + Sbki, the matrix elements Hki undergo a diffusion dynamics along with a finite drift, 
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where Xki-, s = 1 — jgklhki;s with gu = 2 — 5ki and gki = 1 + Ski- The 7 is an arbitrary parameter, giving the variance 
of the matrix elements at the end of the evolution [2]. The above equation describes a multi parametric flow of 
matrix elements from an arbitrary initial condition, say H - However, as discussed in [2], it is possible to define a 
"complexity" parameter Y, a function of various distribution parameters hki- s and bki-, s , in terms of which the matrix 
elements undergo a single parametric diffusion, 
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here fj' implies a product over non-zero 6 W .„ and Xki- S - Further, C is a constant determined by the initial distribution, 
M is the number of all non-zero parameters Xki-, s and bki-, s and f3 = 1,2 for Hamiltonians with and without time- 
reversal, respectively. 

The solution of the eq.(3) gives the state p(H, Y\H , Y a ) of the flow at parameter Y, starting from an initial state 
Hq with Y = Yq. An integration over initial probability density result in the density given by eq.(l) in terms of 
Y(h,b): p(H,Y) = J p(H,Y\H ,Y )p(H ,y a )dH - The evolution reaches a steady state when dp/dY — > with the 
ensemble p{H) approaching the Wigner-Dyson limit, p oc e~(' y W TlH . 

As implied by eq.(3), the variation of p(H\H ), depends on the changes of the parameters hki, bki (for all k, I) only 
through a function Y. This can be proved by considering a transformation of the M non-zero variables of the sets h 
and b to another set {Y,Y 2 , ..,7m} of M variables; hki = hki{Y,Y 2 , ..,1m) and b^i = bki(Y,Y 2 , ..,1m)- As shown in 
[2,5], it is possible to define Y, Y 2 , ..Ym such that the M — 1 variables Y 2l .., Ym remain constant during the evolution 
of p due to any change in sets h, b. The statistics during the transition is therefore governed by Y only. The choice 
of the Y 2 , ..,Ym depends on the system under consideration. For a transition preserving the lattice structure, these 
constants turn out to be the functions of the site-indices in the lattice. For example, the variances hki in the Anderson 
ensemble are functions of the disorder as well as the site-indices k, I; Y can then be identified as a function of disorder 
while Yj (j > 1) as the functions of site indices. Further, as these constants do not appear explicitly in eq.(3), its 
solution and therefore the ensemble-statistics is independent of the specific values of the constants. 

The flow described by eq.(3) can start from any initial state; the only constraint on the choice is that the parameters 
Yj, j > 2, for the initial ensemble should be same as those for the ensemble p(H, h, b). As shown below by an example, 
the initial state can also be chosen as the insulator limit of the disordered system, described by an ensemble of diagonal 
matrices. Although this corresponds to a same value for all initial off-diagonal variances (that is, zero), however a 
choice of different rates of change of hki with respect to Y can result in different possible values for each hki at a later 
stage. 

As an example, consider an Anderson system with a Gaussian site disorder (of variance W 2 /12 and mean zero), 
same for each site, and an isotropic Gaussian hopping with a non-random component (of variance W s 2 /12 and mean 
t s with s = 1,2 for real and imaginary parts respectively) between nearest neighbors (referred as ensemble G later 
on). The corresponding probability density can be described by eq.(l) with 

h kk - W 2 /12, b kk = 

hki-s - fi(kl; s) W 2 /12, b kl , s = f 2 (kl; s) t s 

where fi(kl; s) = 1, f 2 (ki, s) = 1 for for {k, 1} pairs representing hopping, fi(kl; s) — > and f 2 (kl; s) — > for all {k, 1} 
values corresponding to disconnected sites. As obvious, here the distribution parameters hki depend on more than 
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one system parameters, namely, the disorder parameters W, W\ and W2 as well as various functions of site-indices. 
The latter, being invariant of motion, give the parameters Y 2 , ..,Y M . The Y for this case can be obtained by using 
eq.(4), 

N 

Y = a + C (5) 

a = ln|l -jW 2 /12\ 

+ (*/2) £ln [|1 - jW!/6\\t s + 5 t .o\\] + C, (6) 

s 

with M = [3N(N + z(l - 5 m ) + 2 - j3)/2 w (3N 2+e /2. Here zN is the number of connected sites (nearest-neighbors) 
which depends on the topology and the dimensionality d of the system and the e is a function of z, e(z) = (log(iV + 
z(l - 5to) + 2- /3)/k>gJV) - 1; e -» for z << iV. 

Now consider an insulator as the initial state (in the same site-basis as used for G) with zero hopping, that is, 
W s = 0, t s = and a Gaussian site disorder with variance (W 2 /12) = (27) _1 (referred as Go ). This corresponds to 
an ensemble of diagonal matrices with hkk = (27) , hki-, s = for k ^ I and bki- s = for all k,l. A substitution of 
these values in eq.(5) gives the initial value of Y, say Y , where Y — — ^y^o + C with a = — ln2. Note, the basis 
being same, the parameters Yj (for j > 2) are same for both G and Go- (The advantage of choosing the above initial 
state is explained later). As obvious, starting from Go, a variation of diagonal disorder W, hopping parameters W s 
and t s with rates 

§fe = WMkl- »)/6, ^fe = t,h(kl; s) k*l 

Ws 0t 8 

can lead to the ensemble G. Using (dYj/dx) = for j > 2 with x = W, W s ,t s , and, eq.(5) to obtain (dY/dx), it can 
be seen that the above rates correspond to 

^ a 11-7^/121, 
^cxMl-7^/61, 

ay 



cx / 2 t.; (7) 



the variances and means of different matrix elements therefore change with different rates with Y. 

The distribution P of the eigenvalues E n for a metal (for the energy ranges with fully extended eigenfunctions) is 

given by the Wigner-Dyson distribution, P({E n }) = Yli<j \E* — Ej\Pe~* Ek , and, for an insulator by a Poisson 
distribution [7]. The distribution for various transition stages can be obtained by integrating p over the associated 
eigenvector space. Let P({E n }, Y(h, b)) be the joint probability of finding eigenvalues A, of H between Ei and Ei + dEi 
(i = 1, 2, ..,N) at a given h and b, it can then be expressed as P({E n }, Y) = J fliLi — Xi)p(H, Y)AH. Using the 
above definition in eq.(3), it can be shown that the eigenvalues of p(H) undergo a diffusion dynamics along with a 
finite drift due to their mutual repulsion, (see [2] also) 



dP_ x ^ d 

dY ~ ^ ~dE~ n 



dE n E m — E, 



(8) 



Again the steady state of the evolution is given by the limit dP/dY — > 0; P({E n }) in this limit turns out to be a 
Wigner-Dyson distribution. 

The eq.(8) can be used to obtain the correlations between levels. For example, a knowledge of its solution P gives 
the static correlations 
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The P can be obtained by using the analogy of eq.(8) with the equation governing the evolution of the eigenvalues 
of Brownian ensembles (BE) of hcrmitian type [1,3]. The latter, has been studied in great detail in past and many of 
its statistical spectral properties are already known [4]. A brief description of the BE is given in the next section. 

It should be noted here that the single-parametric evolution of the matrix elements of the AE in terms of the 
complexity parameter Y — Yq would result in a similar evolution for their eigenvector components too; this can be 
shown by integrating the eq.(3) over all eigenvalues. However in this paper we confine ourselves to the discussion of 
eigenvalue statistics only; the details for the eigenvector statistics will be published elsewhere. 



The stationary random-matrix ensembles were introduced in the past to model quantum mechanical operators of 
complex systems in which a certain set of quantities (for example, total spin, charge or isotopic spin) was exactly 
conserved; no other integral of the motion existed even approximately [1,3]. The total set of the states of the 
system could then be divided into subsets, each subset corresponding to a particular set of values for the conserved 
quantities. This divides the matrix representation of the operator in various blocks; the deterministic uncertainty due 
to complicated nature of the interactions leads to randomization of the blocks. Due to lack of correlation between 
energy levels of states belonging to different subsets, different blocks are uncorrelated. The statistics of the levels 
within one subset can then be described by a separate random matrix model which can be of various types based on 
the underlying symmetry [1,3]. 

The stationary ensembles arc inappropriate models for systems possessing approximate conservation laws. However 
such systems occur more frequently in practice which motivated Dyson to introduce the Brownian ensembles (BE) of 
random matrices [1,3]. As the latter have been discussed in detail in past e.g. in [3,4] (and references therein), here we 
give only a brief review of the BEs related to Hermitian matrices. Consider the Hamiltonian operator H of a system 
with its elements given by H^i at "time" A and Hki + SHki at "time" A + <5A. A Brownian motion of H is defined by 
requiring that each SH^i is a random variable with the moments: < 5Hki ]S >= — 7-Hfcz ; s£>A, < {SHki- s ) 2 >= guSX [3]. 
The evolution of the distribution of matrix elements, from any arbitrary initial state, can then be given by a Fokkcr- 
Planck equation which has the same form as eq.(3) with 7 a A 2 . For A — > oo, the distribution approaches steady 
state which corresponds to one of the stationary ensembles. The crossover to stationarity is rapid, discontinuous, as 
a function of A, for infinite matrix sizes or very large energy-ranges. 

A Brownian ensemble can therefore be described as a non-stationary state of the matrix elements undergoing 
a cross-over due to a random perturbation of a stationary ensemble by another one. For example, in the case of 
Hermitian operators, a Brownian ensemble H can be given as H = V7(^o + XV) (with / = (1 — A 2 ) -1 ); here V 
is a random perturbation of strength A, taken from a stationary ensemble, and applied to an initial stationary state 
Hq (see also [4]). Using 2nd order perturbation theory, it can be shown that the eigenvalues Ej, j = 1, 2, .., N of H 
execute a Brownian motion too, with their evolution described by an equation same as eq.(8) (with Y oc A 2 /). The 
eigenvalue statistics (e.g. static correlations given by eq.(9)) of a BE can then be obtained by solving eq.(8). The 
eq.(8) is equivalent, under a Wick rotation, to the Schrodinger equation of the Calogero-Sutherland Hamiltonian; 
the equivalence has been used to obtain the eigenvalue correlations for many BEs [4] . It is shown moreover that the 
crossover in correlations is governed, for small A and large N, by a rescaled parameter A which measures locally the 
mean-square symmetry breaking matrix element in units of the mean eigenvalue spacing of H . 

The type of a BE, appearing during the cross-over, depends on the nature of stationary ensembles H 0} V and their 
different pairs may give rise to different BEs [4]. The present knowledge of ten types of stationary ensembles [3] 
leads to possibility of many such cross-overs and, consequently, many types of BEs. For example, the Hamiltonian 
of a disordered system or an autonomous chaotic system, with time-reversal symmetry, can usually be modeled by 
Gaussian orthogonal ensemble. The breaking of time-reversal symmetry e.g. by switching of a magnetic field, with 
A as a measure of the breaking, perturbs the Hamiltonian Ho. The statistical behavior of the system now depends 
on the energy-range of interest. At asymptotically large energies, the statistics can be modeled by Gaussian unitary 
ensembles; however at intermediate energies with sufficiently small values of A, an intermediate statistics (a BE 



II. SPECTRAL PROPERTIES OF BROWNIAN ENSEMBLES 
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between GOE and GUE) would be obtained, indicative of a non-equilibrium behavior. Similarly if the system is 
integrable with regular classical motion for A = and fully chaotic for A ^ 0, the statistics undergoes the Poisson — > 
GOE crossover; (the BE in this case is a superposition of Poisson and GOE ensembles). For many type of crossovers, 
beginning from various stationary states e.g. GOE, GSE, 2GOE, Poisson, uniform etc and approaching GUE in limit 
A — > oo, the 2 nd order correlation functions for all A have been explicitly evaluated [4]; for the other transitions the 
correlations are given implicitly by a hierarchic set of relations [4,2]. 

Here we discuss only the BEs appearing during a transition from Poisson — > Wigncr-Dyson ensemble (referred as 
Wigner-Dyson transition or WDT) caused by a perturbation of the former by the latter (that is, taking Ha, V as 
Poisson and Wigner-Dyson ensemble respectively) . As this transition results in a change of localized eigenstates to 
delocalized ones, its relevance for the study of MIT is intuitively suggested. The BEs related to the Poisson — > Wigncr- 
Dyson transition can be described by a N x N ensemble H represented by the following probability distribution for 
all (independent) matrix elements: 



p(H) oc exp 



N 

<]>>!-2 7 (i+ M )]ri!^i 2 

i—l i<j 



(10) 



with (1 + fi) = (A 2 /) -1 ; here H = H for A — > or /x — > oo. An ensemble H given by the above measure, is also known 
as Rosenzweig-Porter ensemble (RPE) ; Note it also corresponds to an ensemble of Anderson Hamiltonians with very 
long range, isotropic, random hopping. 

The eq. (8) describes the evolution of the eigenvalues of a generalized Gaussian ensemble with a probability density 
(1) and is therefore applicable for the BEs defined by a probability density (9) too. A comparison of measure (9) 
with measure (1), gives a variance hki- s = (47(1 + m)) -1 ; h^k-s = (2 7 ) _1 and mean bki- s = for all fc, I and s indices. 
Using these values in eq.(4), the parameter Y for the BE case can be given as 

~-^ + Fo (for /i»l) (11) 
^7M 

with M = f3N(N + 2 — j3)/2 and Y a = ^pj m 2 + C as the complexity parameter of the ensemble H (note, Y = Y a 
for [i — ► 00). 

A typical matrix in the ensemble (9) has the diagonal elements of order j^ 1 / 2 and off-diagonals of the order of 
( 7 /i) _1 / 2 (= o(Y — Yq) 1 / 2 ). The number of off-diagonals being N times more than the diagonals, the matrix behavior 
is governed by the parameter /x. Thus, for large BE (N — > 00), a radical change from Wigner-Dyson case can only 
occur if n increases more rapidly than N (which makes the total strength of the off-diagonals weaker than that of 
diagonals) . This results in three different regimes of the mean-level density R\ (E) [9] : 

N ^2 

R 1 (E) = ^c- E for N(Y - Y ) -» (12) 

\ 7T 



= NF(E, a) for N(Y -Y )=a (14) 

with a as an A-indcpcndcnt constant. Although the exact form of the function F{E, a) is not known, its limiting 
behavior can be given as follows: F(E,a) w eT^ /y/n for a « 1 and F(E,a) w (47r7a) _ 1 \J 87a — E 2 for a >> 1, 
E 2 «a[9]. 

The 1 st order correlation Ri, also known as mean level density, changes from an exponential to semi-circular form 
at the scale of (Y — Y ) ~ NAf with A/ as the local mean level spacing; the evolution of R\ can therefore be described 
in terms of the parameter (Y — Ya). However the transition of higher order correlations R n (n > 1) occurs at a scales 
determined by (Y — Ya) ~ Af [4,9,2]. As a result, their transition to equilibrium, with \Y — Yq\ as the evolution 
parameter, is rapid, discontinuous for infinite dimensions of matrices [1]. But for small- Y and large N, a smooth 
crossover can be seen in terms of a rescaled parameter A(E): 
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A(E,Y) = \Y-Y \/A* 



(15) 



For finite N, A varies smoothly with changing \i: A = R\j^[i. This results in a continuous family of BEs, parameter- 
ized by A, existing between Poisson and Wigner-Dyson limit. However the level-statistics for the large BE (N — > oo) 
can be divided into three regions [9]: 

(i) Poisson regime: N 2 (Y — Yq) — > 0: The off-diagonals, responsible for the correlation between levels, are 
negligible. The lack of repulsion between levels results in a mean level spacing A; oc N~ x (see eq.(12)), thereby, giving 
A — > and the Poisson statistics. 

(ii) WD regime: N 2 (Y — Y ) — > oo: The contribution from both, the diagonals as well as off-diagonals is of the 
same order, leading to long-range correlations between levels. The repulsion of levels now results in a mean level 
spacing A; oc N~ x / 2 . (see eq.(13)) which gives A — > oo and Wigner-Dyson statistics. 

(iii) Critical regime: N 2 (Y — Yq) = (A'yc)^ 1 — a constant: For [i = cN 2 with c as a constant independent of 
N, a sequence of approximately o{l/y/c) levels show Wigner-Dyson behavior. The more distant levels display weak 
correlations of the type existing near the Poisson limit resulting in a A; k o(\/N). The parameter A is therefore 
Af- independent: 



note it is also independent of the symmetry parameter /3. 

The finite, non-zero A-valuc for fj, — cN 2 in limit N ^ oo therefore gives rise to a third statistics, intermediate 
between Poisson and Wigner-Dyson ensemble, which is known as the critical Brownian ensemble (CBE) . This being the 
case for arbitrary values of c (non-zero and finite), an infinite family of critical BE, characterized by c (or \i c = cN 2 ), 
occurs during WDT. Note that the critical BEs, with c — > oo and c — > 0, correspond to the Poisson and Wigner-Dyson 
limit, respectively. 

The presence of a family of the critical BEs can be seen from any of the fluctuation measures for WDT. One 
traditionally used measure in this regard is the relative behavior of the tail of nearest-neighbor spacing distribution 
P(s,A), defined as a(d,A) — J Q 5 (P(s, A) — P w (s))ds/ J^(P p (s) — P w (s))ds with 5 as any one of the crossing points 
of P w (s) and P p (s) (here subscript w and p refer to the Wigner-Dyson case and Poisson case respectively) [10]. In 
the limit N — > oo, a — and 1 for Wigner-Dyson and Poisson limit respectively. The figure 1 shows the numerically 
obtained behavior of a (for S w 2.02) with respect to \z — c\ (=\fj, — n c \N~ 2 ) for a fixed c (arbitrarily chosen) with 
z as a variable; Here z and c are the values of the parameter /iN~ 2 for a general BE and a critical BE respectively. 
The constant value of a at \z — c\ = for different N- values confirms the size- independence of the level-statistics of 
BE with parameter [i = cN 2 and therefore its critical nature; we have verified it for other c-values too and find, for 
finite, non-zero c-values, < a < 1. Further the convergence of a- values for BEs with different fi and N- values on 
two branches indicates the presence of a scaling behavior in the level-statistics of BEs with \z — c\ (=\(i — ii c \N~ 2 ) as 
the scaling parameter. 

As shown in figure 1(a), a for a critical BE is between and 1. A fractional value of a indicates a tail-behavior 
of critical BE different from that of Poisson as well as Wigner-Dyson limit. As shown in figure 1(b), the P(s) for a 
critical BE with a finite parameter c has an exponential tail, P(large s) ~ exp(— us); this behavior of P(s) is also 
sometimes referred as semi-Poisson distribution, due to presence of repulsion at small energy-scales and exponential 
decay at large separations. 

III. ANALOGY BETWEEN BROWNIAN ENSEMBLES AND ANDERSON ENSEMBLES 

The same evolution equations of P for AE and BE imply a similarity in their eigenvalue distributions for all Y- values, 
under similar initial conditions (that is, P(/x, Y ) same for both the cases). As a result, one obtains the analogous 
evolution equations for their correlations R n too. The mean level density R\(E,Y — Yq) of an AE can therefore 
be given by the level-density of a BE with same \Y — Yq\ value (and appearing during a Poisson to Wigner-Dyson 
transition). Similarly, the analogy of evolutions of R n (n > 1) in the two cases implies 

(i) a smooth crossover of R n for finite size Anderson systems in arbitrary dimensions, 



A(E) = (l/4c7T7)e 



(16) 
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(ii) the parameter A governing the smooth crossover of R n for finite size AEs can again be defined by eq.(f5), with 

Y — Yq given by eq.(4) and A; as the local mean level spacing for AE (see [11] also), 

(hi) the correlations R n , n > 1, of an AE can be given by those for a BE with a same A value although their 
parameters Y (as well as level densities) may be different; 

(iv) the discontinuity of the transition of R n for infinite size of Anderson matrix, 

(v) the existence of a size-invariant level statistics, different from two end-points, if an AE has A = 
size — independent; the statistics survives the thermodynamic limit L — > oo. As explained later by an example, 
the above condition on A is satisfied at the critical point of d > 2 dimensional Anderson Systems; the corresponding 
level-statistics is referred as critical. . 

The implications (i,iv) are well in agreement with known results about AE-correlations [13]. The implications 
(ii,iii,v) indicate the single parametric dependence of the level statistics for AEs. The parameter A for the AE and 
BE will henceforth be referred as A a and A&, respectively. The level-statistics of a finite-size AE at A a is then given 
by a BE with its parameter \i satisfying the condition 

A a = A b (17) 

where A = Rf/A-y/j, with i?i as the level-density of the BE; the determination of A a is explained later by using an 
example. As the BEs with different combinations of the parameters /i and N can have same Ab, the correlations of a 
finite size AE can be mapped to many BEs. However the critical BE corresponding to a critical AE is unique; this 
can be understood as follows. Using eq.(16) in eq.(17), the parameter c for a critical BE corresponding to an AE can 
be given by 

c= (47r 7 A a )- 1 e- B2 . (18) 

The A a for a critical AE being size- independent, its critical BE analog remains same for all system sizes. However, 
the A a for an AE, away from its critical point, is size-dependent and therefore corresponds to different c values (that 
is, different critical BEs) for different system sizes. 

The A a for a disordered system can be determined by a knowledge of Y — Yo and A; . The complexity parameter 

Y — Yq is system-specific and depends on various system parameters. For a <i-dimensional disordered system of linear 
size L, the local mean level spacing A; within a correlation volume of linear dimension £ is related to mean level 
density R\. A; = (L/^R^ 1 where £ is the localization length or correlation length in case of localized states and 
extended states, respectively [12,13]. The £ can be determined by a knowledge of the wave-function correlations e.g. 
inverse participation ratio I2 [17]: ( d oc (/2) -1 for localized eigenstates [17]. As mentioned in the last paragraph of 
section I, the wave-function statistics and, therefore £, can also be described, in principle, by a complexity parameter 
formulation. However, the related work being still in progress, we use, in this paper, the £ results given by previous 
studies. 

Let us consider the example G given in section 2; its parameter Y is given by eq.(5). The initial state Go has 
the parameter Yo same as that of the initial state chosen in the BE case in Section II. Note as N\Y — Yo\ - * an 
A^-independent function for the case G, its Ri is given by eq.(14) with Ri = NF. Using eq.(5) in eq.(15), A for the 
case G can be given as 

A a (^,y) = (^^)c 2 ^- d (19) 

with F(E) giving the energy-dependence of A (as |e| w for large N). Following eq.(17), the level-statistics at 
A a — > 0,oo corresponds to Poisson (or insulator limit) and Wigner-Dyson behavior (metallic limit), respectively. In 
finite systems, a change of disorder results in a smooth variation of £ as well as a — ao and, therefore, A a which 
induces a crossover of the level- statistics from Poisson — > Wigner-Dyson ensemble. The intermediate level-statistics 
at each A a of a finite-size AE is then given by a N\ x N\ BE with its parameter /j, satisfying the relation A a = A : 

Ii ki l3{Air\a-a \F 2 )- l C 2d L d Rl (20) 

with Ri = Ri(E; n, N\) as the level-density of the BE. As obvious, the determination of n from the above equation 
is not easy, its both sides being ^-dependent. The Ri for the critical BEs being ^-independent (given by eq.(12)), 
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it is preferable to map an AE to a critical BE; the substitution of eq.(19) in eq.(18) gives the c-parameter for the 
corresponding critical BE: 

c«/3(47r|a-ao|^ 2 c £;2 )" 1 r 2d i d (21) 

Thus each state of disorder in an AE of size L (N — L d ) can be mapped to a critical BE with the parameter c given by 
eq.(21). Note, the right side of eq.(21) being energy-dependent, different energy ranges of a given AE will, in general, 
correspond to different critical BEs. 

Equation (19) indicates the sensitivity of the parameter A to localization length Q and system size L. It is now well- 
known that ( is a function of disorder-strength, energy, system-size L as well as the dimensionality of the system. For 
systems with finite L (in arbitrary dimensions d > 1), at a fixed energy, decreases with increasing disorder-strength. 
Consequently, in the strong disorder limit (where £ ~ o(L° j), A a — > and the level-statistics of the AE approaches 
Poisson behavior (as A& — > for its BE analog). In the opposite limit ( <~ o(L) of weak disorder, A a — > oo and, 
therefore, the statistics of the AE is given by a BE at A& — > oo which corresponds to Wigner-Dyson behavior. By a 
suitable choice of disorder, however, it is possible to achieve finite values of the ratio ( 2 /L (due to finite L) in arbitrary 
dimensions which in turn gives finite, non-zero A a and, thus, a finite c for its BE analog. The latter implies that the 
AE-statistics is intermediate between Poisson and Wigner-Dyson limit, with an exponential decay of the tail of its 
nearest-neighbor spacing distribution P{s). For finite L, therefore, a smooth crossover from Poisson to Wigner-Dyson 
statistics can be seen, for any dimensionality d > 1, as a function of A a by varying the disorder-strength. Note, two 
finite size AEs of different dimensions can show same level-statistics if their parameters A a are equal. For example, 
consider the behavior of levels of a one-dimensional AE of size L and at a disorder strength which gives Y = Y\ . The 
behavior will be same as that of a three-dimensional AE of linear size L, at a disorder strength which gives Y = Y3 
where Y 3 = Y^R^ /(sR^) 2 ; here R[ d) and Q refer to the mean level density and localization length in dimension 
d, respectively. 

The dimensionality dependence of the Anderson transition and the critical level-statistics is well known. For 
example, the level statistics at the critical disorder for d > 2 dimensional, finite systems shows a " semi- Poisson" 
behavior which survives the infinite size limit. The same behavior is seen for d < 2 dimensional finite systems, in a 
regime where ( <~ L, however the statistics approaches a Poisson behavior in the thermodynamic limit. The above 
behavior can be explained within " A-formulation" . As mentioned in the section II, a " semi- Poisson" behavior of 
the level-statistics is a characteristic of critical BEs with finite c parameters and therefore of the AEs with finite A a 
parameters (see eq.(19)). The AE-statistics is expected to maintain its semi-Poisson behavior even in thermodynamic 
limit if A a =size- independent. In this sense, A a can be identified with the dimensionless conductance g\ both g 1 A a ► 
0, 00, constant correspond to same statistical limit, namely, Poisson, Wigner-Dyson and a critical level statistics, 
respectively. In fact, A a can be expressed in terms of the dimensionless conductance g of the system. This is because 
g is connected to £ (based on scaling theory of localization for disordered systems [14]): £ oc Llogg -1 for exponentially 
localized states, ( oc \{g/g c ) — 1|~" near the critical point and ( ~ o(L) with g oc L d ~ 2 . For example, using eq.(19), 
the A a — g relation for the case G near the critical point can be given as: 

K{E,Y) = ( '"'y 2 ^ ) \{g/g c ) l\-™L~ d (22) 

with g c as the critical point conductance and v as the critical exponent. 

As indicated by eq.(19), the size-independence of A is governed by the size-dependence of the localization length. 
For example, for a d-dimensional disordered system, with \Y — Yo\ ~ o(L Xld ), the level-density is Ri(Y — Y ) ~ o(L X2d ) 
where 

x 2 = 1 for x\ < —1, 

x 2 = (1 - x{)/2 for x x > -1. (23) 

Thus the critical point of the level-statistics (that is, A=size-independcnt) can exist only if, in thermodynamic limit, 
the disorder conditions in the system give rise to a localization length ( ~ o(L Xs ) where X3 w (2 — x\ — 2x 2 )/2 or, 
equivalently, 
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x 3 = kil/2 for X\ < — 1 

x 3 = 1/2 for xi > -1. 



(24) 



The existence or non-existence of a critical level statistics in an AE therefore depends on the the size dependence of 
( which in turn is sensitive to the dimensionality of the system [14]: 

Case d < 2: For a d = 1 disordered lattice, almost all states are known to be exponentially localized even in a weak 
disorder limit. The ( in this case is finite, C ~ nl <~ o(L°) [14] (with I as the mean free path), which gives x 3 = 0. 
As obvious, the condition (21) can not be satisfied for any x 1; equivalently, for any Y — Y (e.g. for any disorder 
conditions). As a consequence, a critical level-statistics can not occur in one-dimensional case. 

Equation (16) suggests that a " semi-Poisson" type statistics can be seen in d = 1 case for L of the order of few 
mean free paths (i.c for finite l 2 /L). In limit L — > oo, however, A a — > and the level statistics approaches Poisson 
behavior irrespective of the disorder strength. 

In two dimensions, the perturbative estimate of the localization length is £ rs ^exp[-7rfci?Z/2] with kp as the Fermi 
wave number [14] and, in the limit L — > oo, electronic states are expected to be localized even for small microscopic 
disorder [15]. This again corresponds to x 3 — (thus absence of critical level-statistics), and, the Poisson statistics 
for the levels in thermodynamic limit (as A a — > 0). Note, however, that due to exponential nature of £, the ratio 
( 2 /L can be kept non-zero and finite (by changing disorder) for a large range of L. The system can therefore show 
the semi-Poisson statistics in a large range of system sizes. 

Case d > 2: For d > 2 dimensional, infinite systems, the change of disorder W leads to a discontinuous change in 
£ and thereby A a : ( cx |1 — (W/W c )\~" . Here v is the critical exponent and W c is the critical disorder. For W > W c , 
almost all states are exponentially localized with ( <~ o(L°) which results in A a — > and Poisson behavior of the 
statistics. For W < W c , the derealization of states occurs with ( — > o(L); this gives A a — > oo and the Wigner-Dyson 
statistics. At W c , however, the inverse participation ratio I 2 for d > 2 case shows an anomalous scaling with L [16]: 
I 2 oc L~ D2 with D2 as the multifractality exponent. This gives £ d oc< I 2 >~ 1= (qL D2 or £3 = 1)2/^ with Co as a 
size-independent function; note D 2 = D 2 at the critical point [18,19]. The size-independence of the level statistics at 
the critical disorder therefore requires 

D 2 = d\ Xl \/2 for X!<-1 

= d/2 for X!>-1. (25) 

For example, as x\ w —1 in the case G, the existence of its critical point requires D 2 ~ d/2. Note the numerical 
results for D 2 , at the critical point of a d = 3 dimensional AE system (of type G) fluctuate in the range 1.4-1.6 
[20,22-24] (also see references [77,79] in [23]); this is in close agreement with the result given by cq.(25) for d — 3 case. 
The above prediction for D 2 can be used to determine the critical BE analog for the critical state of the AE example 
G for d > 2 case: 

cw (4n\a-a \(: 2 d F 2 e E2 )- 1 p (26) 

(as e w 0). Thus, unlike the d = 1 case showing only Poisson level-statistics in the thermodynamic limit, the energy 
levels of an infinite size AE for d > 2 case can show three types of behavior, namely, Poisson, Wigner-Dyson and a 
critical BE type statistics, at the disorder strength above, below and at the critical disorder, respectively. 
The study [22] suggest a connection between D 2 and the level-compressibility \: 

D 2 = d(l - 2 X ) (27) 

A comparison of eq.(25) with eq.(27) gives the \ f° r d > 2 dimensional AE at the critical point: \ w 0.25. (Note the 
above x-result is valid only for the cases of type G with ( oc L D2 ^ d ). The tail of the distribution P(s) is also believed 
to be related to D 2 [26]: P(large s) rts e~ KS where k — {2\)^ 1 ~ 2. The results for \, k are in close agreement with 
earlier numerical studies on Anderson systems [20,22-25]; our numerical study, given in section V, also confirms the 
above results. The symmetry-independence of our theoretical prediction for \ and n for Anderson systems is also 
in agreement with numerical observations [25,24]. As discussed later, however, the eq.(27) (and therefore above x> k 
results) seems to be valid only in the weak multifractality limit i.e D 2 <~ d (see paragraph below eq.(36)). 
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For disordered systems, in general, both Y — Y as well as £ are a function of coordination number, disorder strength, 
hopping rate, dimensionality as well as boundary conditions of the lattice. The changing complexity due to change 
of the system parameters plays the role of a random perturbation, of strength \JY — Y , applied to the system. Here, 
again, the statistics of the levels is governed by A and, therefore, by the competition between local mean-level spacing 
A; and the perturbation strength Y — Yq. The perturbation mixes fewer levels with increasing system size if A; 
increases with L at a rate faster then that of yF — Yq and, as a consequence, leaves the level-statistics unperturbed 
in limit L — > oo. In the opposite case with slower rate of change of A; with L (as compared to \JY — Yq), even a 
small change in the complexity parameter is capable of mixing the levels in an increasingly large energy range of many 
local mean level-spacings. This results in an increasing degree of the eigenfunctions derealization and Wigner-Dyson 
behavior of level-statistics in thermodynamic limit. The critical regime occurs when both \JY — Yq and A; change 
at a same rate with L; the perturbation in this case mixes only a finite (non-zero), fixed number of levels even when 
the system is growing in size. As A remains finite in limit L — > oo, it gives rise to a new statistics different from 
the two end-points (A — > and oo). Note, the disordered systems with different dimensionality can have different 
critical values of A (due to dimensionality dependence of A; as well as \Y — Yq\) and, therefore, correspond to critical 
BE analogs with different c values. Further the boundary conditions/ topologies, leading to different sparsity and 
coordination numbers, can also result in different critical level statistics even if the underlying symmetry and the 
dimensionality is same; this is in agreement with numerical observations [20] and analytical study for 2D systems [21]. 
A knowledge of A can then be used to map the critical level statistics at MIT for various dimensions d > 2 — > oo to 
the infinite family of critical BEs. 

IV. DETERMINATION OF FLUCTUATION MEASURES FOR MIT 



Many results for the spectral fluctuations of the WDT with Poisson ensemble as an initial state are already known 
[4] and can directly be used for the corresponding measures for the MIT in different disordered systems. 

A. MIT With No Time-Reversal Symmetry 

The fluctuation measures, for the Anderson transition in presence of a magnetic field, can be given by the BEs 
appearing during a WDT which violates time-reversal symmetry. Such a WDT, occurring in a complex Hermitian 
matrix space (that is = 2), corresponds to a transition from Poisson — ► GUE ensembles. 

(i) The 2-Level Density Correlator R 2 (r;A): The R 2 for BEs during Poisson — ► GUE transition has been 

obtained by various studies [4,27,28]. Here we use the form given in [27] for the purposes to be explained later (note, 
our A is equivalent to A 2 /2 used in the [27]), 

A A /*0O 

R 2 (r; A) = 1 + — / du F c -^ 2 -4.Au (28) 
r Jo 

with 

F = sin(w)/i — cos(ur)/ 2 

h = {2/z)[h{z)-x/^h{z)} 

h = (l/u)[I 2 (z) - y/2^h(z)} (29) 

where z = \/32irA 2 u 3 and /„ as the n th Bessel function. (Note, the eq.(4.15) in [27] has a misprint in the coefficient 
of u is the exponent; the correct coefficient is given in the eq.(28) above ). 

The eq.(28) gives the exact form of two-point correlation for the Anderson transition with no time-reversal symmetry. 
Here R 2 (r, oo) = 1 — (sin 2 (7rr)/-7r 2 r 2 ) and R 2 (r, 0) = 1 corresponding to metal and insulator regime respectively. A 
substitution of critical value of A a in eq.(28) will thus give R 2 for the critical AE. 
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For large A- values (for all r), R 2 can be approximated as follows [27,28]: R 2 = 1 — Y 2 where 

-4A 1 r 2 

Mr, A) = 16 ^2 A2 + r2 ~ ^v^Me" 5 * - 1] (3°) 



3 sin (nr) 



2^ 2 A s i n h 2 (rV372A) 



(for r « VA) 



However, for r > VA, Y 2 (r,A) = - 167r ^ +t ^ + 2iT l ri . As A = (Acirj) 1 (near E = 0) for a critical BE, the 
Y 2 ss (1 - 87r 2 A)/27r 2 r 2 for r > 2/?ttA (here = 2). 

The above large r-behavior of l^O"; A) at A = A;, results in a non-zero, fractional value of the sum 1= Y 2 (r; A)dr 
for a critical BE of complcx-Hermitian type: 

I ml -(fin 2 A)- 1 . (31) 

Note a < I < 1 value is believed to be an indicator of the multifractality of the wavefunctions and the fractional 
compressibility of the spectrum; (J = 1,0 for the WD and Poisson case, respectively) [22,29]. A fractional behavior 
of I and the multifractality is already known to exist in critical AE [22,29]. Using A a = A;, in eq.(31), one can now 
determine the measure I for an AE: J w 1 — ((3w 2 A a )~ 1 . 

(ii) Nearest Neighbor Spacing Distribution P(S) 

The nearest-neighbor spacing distribution P(s) for the MIT with no time-reversal symmetry can similarly be given 
by using the one for the BE during Poisson — > GUE transition [30,31]: 



P(s; A) oc _^ c - 2 /8A f°° dx e^M ^'W 
V2^A Jo x 



(32) 



A substitution of A — > oo and A — > in the above equation gives the correct asymptotic limits, namely, Wigner-Dyson 
and Poisson, respectively: P(s; oo) = P w (s) = 32s 2 c~ 4s / w /tt 2 (WD limit) and P(s;0) = P p (s) oc c~ s (Poisson limit). 
(Although this result is rigorous for 2 x 2 matrix space but is proved reliable for systems with many levels; see [31]). 

(iii) Level-Compressibility 

The level compressibility \ = 1 — Y 2 (r)dr = 1 — / is an important characteristic of the critical level statistics 
and the multifractal nature of the wavefunctions. 
The x f° r a BE can be obtained by using eq.(28), 



/•OO 

X(A) = 1 - 4ttA / du/!(z)exp[-2Au 2 - 4ttAu] 
Jo 



(33) 



w 1 - 4tt 2 A for small A (34) 
w (2tt 2 A)- 1 for large A (35) 

The substitution of eq.(16) for At in eqs. (29,30) gives, in the band around E = 0, \ = 1 — (n/jc) and \ rts (A^c/2-k), 
in small and large A limits, respectively. Thus a critical BE characterized by a finite c value shows a fractional 
level-compressibility. As clear from the above, \ — * 1 for A — > (or c — > oo) which corresponds to a Poisson behavior, 
and, x ^ for c ^ or A ^ oo which corresponds to the GUE statistics. 

The compressibility of the energy levels of Anderson systems at their critical point is already known to be fractional, 
with x = 0, 1 in the metallic and the insulator phase, respectively. The existence of a fractional x f° r both critical 
BE and critical AE is consistent with our claim about their spectral analogy. The compressibility of the AE with 
different types of disorders and lattices can now be obtained just by finding the same for their critical BE analogs. 

For the critical BE case (d — 1) with large Afc (cquivalently, small c), eq.(27) along with eq.(35) gives 

D 2 = 1 — Ajc/ir for small c (36) 

The eq.(36) gives the correct fractal dimension in the limit c — > (the Wigner-Dyson limit): D 2 (c = 0) = 1. However, 
for small A& (or large c ), eq.(27) implies D 2 = and therefore D 2 = — 1 in the Poisson limit c — > oo, which is 
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different from the expected result D 2 = for the localized states. As mentioned in [18], similar violation of eq.(27) is 
indicated by numerical data for the tight- binding models in dimensions d > 4. The observed inaccuracy of eq.(27), for 
both AE as well as BE in the strong multifractality limit, also lends credence to our claim regarding AE-BE analogy. 

It is worth mentioning here that, similar to the AE-BE mapping, the spectral statistics of any generalized Gaussian 
ensemble, with probability density given by (1), can be mapped to BEs [2]. The non-validity of eq.(27) for critical 
BEs with small c parameters implies, therefore, the same violation for all generalized Gaussian ensembles in strong 
multifractality limit. The implication is already known to be correct for the power law random banded matrices 
[18,19,35] (also see section VI) and for the random matrix ensemble introduced by Moshe, Neuberger and Shapiro 
(later referred as MNS model) [34]. 



B. MIT With Time- Reversal Symmetry 

The statistical measures for the Anderson transition in presence of a time-reversal symmetry can similarly be 
obtained by using their equivalence to a WDT preserving the same symmetry, that is, a transition from Poisson 
— > GOE ensembles; the later occurs in a real-symmetric matrix space (here (3—1). However due to the technical 
difficulties [4]), only some approximate results are known for the latter case. 

(i) The 2-Level Density Correlator R 2 (r;A): 

The R 2 for small-r can be obtained by solving eq.(17) of [2] for = 1 which gives 

R 2 {r,A) w (7r/8A) 1/2 re- r2/16A / (r 2 /i6A) (37) 

with Jo as the Bessel function. 

Similarly for large-r behavior, R 2 can be shown to satisfy the relation (see eq.(23) of [4]) 

R 2 (r,A) w R 2 (r, oo) + 

J-oo (S 2 + 4tT 2 /3 2 A 2 ) 

« R 2 (r, oo) + 20A/(r 2 + 4tt 2 /? 2 A 2 ) (39) 

where 0=1 and R 2 (r, oo) = 1 — sin 2 (7rr)/-7r 2 r 2 — (J r °° da;sin7nr /Va;) (^sin7rr/7rr) (GOE limit). 

As can be seen from the above, Y 2 (r, A) w — 4 ^i^i +r i + 27T 2 r -i for r > \/A. However, note, for r > 27rA, the behavior 
of Y 2 is different from that of a BE with no TRS: Y 2 w (1 — 4/37r 2 A)/27r 2 r 2 . This further suggest following behavior 
of I: 1=1 — (0tt 2 A)~ 1 . The I for a critical BE is therefore symmetry-dependent; (as A = A& docs not depend on 0). 
However the I for its AE analog is independent of the symmetry parameter 0; this is because A = A a oc 0^ 1 in this 
case (see eq.(19)). 

(ii) Nearest Neighbor Spacing Distribution P(S) 

The P(s) for this case can be given by using the one for a BE during Poisson — >GOE transition [30]: 

P(s,A) = (V8A) 1 / 2 se- s2/16A / (s 2 /16A) (40) 

with Jo as the Bessel function; note, as expected, this is same as R 2 behavior for small-r (eq.(37)). 

(iii) Level-Compressibility 

The lack of the knowledge of R 2 (r, A) for entire energy-range handicaps us in providing an exact form of the 
compressibility for the time-reversal case. However its approximate behavior can be obtained by using the relation 
X = 1 — I- Thus, for a time-reversal critical BE (0 = 1), 

X «(tt 2 A)- 1 (41) 
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and therefore \ w (7r 2 A a ) _1 for its AE analog. 

Equations (35, 41) indicate the influence of underlying symmetry on the compressibility of the levels: x ~ (/37r 2 A) _1 . 
Note % for a BE is symmetry-dependent due to A = Aj being /3-independent (see eq.(16)). However as A = A a cx 
for an AE (see eq.(19)), its \ would be symmetry-independent; this is in agreement with numerical observations 
for Anderson systems [25,24,36]. This further implies that the critical BEs corresponding to critical AEs with and 
without time-reversal symmetry would be different. 

In past, an attempt to explain the symmetry independence of the level-statistics at the Anderson transition was 
made in the study [36] by suggesting a scaling behavior of the distribution P(s) with the conductance g and the 
symmetry parameter [3. The P(s) in the study [36] was obtained by interpolation between metallic and Insulator 
limits. We note that, by using the A — g connection (eq.(22)), the P(s) given by eqs. (32,40) can also be expressed as 
a function of g. 

V. NUMERICAL COMPARISON OF THE LEVEL-STATISTICS OF CRITICAL AE AND CRITICAL BE 

In this section, we investigate the AE-BE spectral analogy by numerically comparing two of their fluctuation 
measures, namely, P(s) and the number variance £ 2 (r) =< (r— < r >) 2 >. The former is a measure of the short- 
range correlations in the spectrum and the latter, describing the variance in the number of levels in an interval of r 
mean level spacings, contains the information about the long-range correlations [3]. The £2 is also an indicator of the 
compressibility of the spectrum; lim r — > 00 £2 (r) ~ X r - To study the AE-BE analogy in presence of time-reversal 
symmetry as well as its absence, we consider two cases of the 3-dimensional AE (simple cubic lattice of size L = 13 
and with Gaussian site disorder) in critical regime: 

(i) AE t : 

The AE with isotropic random hopping, hard wall boundary conditions and time-reversal symmetry; here W — 4.05, 
W\ = 1, W2 = 0, t\ = 0, ti = 0. The criticality of AE for same disorder parameter values but with periodic boundary 
conditions is numerically studied in [32] . However the system remains in the critical regime under hard wall boundary 
conditions too. 

(ii) AE„t : The system G with isotropic non-random hopping, periodic boundary conditions and no time-reversal 
symmetry; here W = 21.3, W\ = 0, W2 = 0, t\ = 1, ti = 1. The time-reversal symmetry is broken by applying 
an Aharnov Bohm flux (f> which gives rise to a nearest neighbor hopping Hki = exp(i(f>) for all k, I values related to 
the nearest-neighbor pairs [24]. The flux <p is chosen to be non-random in nature, that is, < cos 2 (</>) >= W\ = 0, 
< sin 2 (</>) >= W 2 — and < cos(^) >= h = 1, < sin(0) >= t 2 = 1. 

We study each AE case for two system sizes L = 10 and L = 13 by numerically diagonalizing the matrices of 
the ensembles by standard techniques. Each ensemble contains few thousand matrices and the statistical average is 
performed approximately over 3 x 10 5 levels, obtained by taking 200 levels in a small energy range around center 
E = of the spectrum of each matrix. Each BE (chosen with 7 = 2) is also analyzed for two dimensions N = 1000 
and 3000. Note, due to a higher rate of change of the mean level density, the A (eq.(16)) for BEs changes more rapidly 
with energy as compared to AE cases (eq.(19)). To avoid mixing of levels with different transition rates, therefore, 
fewer (w 100) levels are taken from the spectrum of a matrix in the BE case; the total number of levels for BE 
analysis is kept nearly same (as in the AE cases) by taking a bigger ensemble. Each spectrum is unfolded for P(s) 
and S 2 (^)-analysis. The unfolding is carried out by numerically calculating the unfolded levels rj = Jjf^ dxRi(x) for 
(j = 1,2,.., N) with symbol Ej used for levels before unfolding. 

The parameter Y — Yq for both AE cases is given by eq.(5). As N\Y — Yq\ ~ o(l), the mean level-density R\ for 
both AE cases is given by eq.(14); R\ — NF with F as an N- independent function of energy. This is confirmed by 
our numerical analysis of Ri for different N- values for each AE case, with figure 2(a) showing the comparison only for 
two A^-values. (Note, for R\ study, the spectrum is not unfolded and almost all eigenvalues of each matrix are used 
for the analysis). As mentioned in section III, the critical BE analogs for the fluctuations measures need not have a 
same R\. The i?i-behavior for the critical BE analogs for the fluctuations of the AE cases is shown in figure 2(b); 
the numerical fitting confirms that R\ = (N / ^/tt)c^ e for each critical BE case which is quite different from their AE 
analogs. 

The figures 3 and 4 show the P{S) and T,2(r)/r behavior for the two AE cases. The almost same behavior for two 
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system sizes in each AE case, for both the measures, confirms their critical nature. We find, from figure 4, that the 
large-r behavior of E 2 (r)/r for both AEs seems to converge to \ ~ 0-25 which confirms our analytical prediction of 
X for AEs (with ( ~ o(L D ' 1 l d ) 1 based on equality of eqs.(25,27)); it is also in agreement with other numerical studies 
[20,22-24]. 

The determination of the critical BE analog of the fluctuation measures of an AE requires a prior information about 
A a (given by eq.(19)). Although we know the function F for each AE case (see figure 2(a)) as well as a — ao (from 
eq.(6), a — ao = 1.36,5.43 for AE t and AE nt , respectively), however the determination of £ requires an statistical 
analysis of wavefunctions. Fortunately equations (35,41) suggest that the parameter c of the critical BE analog of a 
critical AE can also be obtained (approximately) from its \ behavior: c w (/?7rx/47). Using \ ~ 0.25, the theoretically 
expected c parameters for AE t and AE nt are 0.1 and 0.2, respectively. The figures 3 and 4 confirm the existence of 
the critical BE analogs, of the AE cases, at the above c values. Note the above relation between the parameter c of 
a critical BE and x °f a critical AE is obtained by combining the theoretical results for (i) critical point D 2 behavior 
predicted by our A formulation, (ii) D 2 given by eq.(27), (iii) x f° r a BE (eqs. (35,41)), (iv) AE-BE analogy. The good 
AE-BE agreement observed in figure 3, 4 therefore indicates the validity of all the above formulations, used to derive 
c{BE) - x( AE ) relation. 

The exponential decay of the tail of P(s) for d > 2 AE system at the critical point has been confirmed by various 
numerical studies (for example, see [7,20,22-24,33]). The validity of AE-BE analogy requires a similar decay of P(s) 
for the critical BEs too. The figure 1(b) compares P(s) behavior for a few critical BEs with the function e~ KS . The 
fitted k values are close toK« (/?7r/87c) for intermediate c ranges; the k — c relation is obtained by using k = (2%) _1 
and x — (47c//?7r) at E — (see below eqs. (35)). The insets in figure 3 compare the tails of the P(s) for the AEs 
and their BE analogs with function e~ KS ; we find (swl.5 - 1.7. The result is close to our analytical prediction k w 2 
for the critical AE case (see below eq.(27)). However the lack of exact agreement seems to suggest the approximate 
nature of the x~ K relation, namely, k = (1/2%); (note as our analytical prediction x = ^l c l fl n is found to be in 
excellent agreement with numerics, this leaves x — K relation as the possible source of error). 

. The study [29] claims that the critical level-statistics in the Rosenzweig-Porter ensemble (similar to BE, as 
mentioned in section II) does not have a fractional compressibility and, therefore, is different in nature from that for 
critical AE. However our analytical results, supported by the numerical evidence, disprove their claim. Our numerical 
study confirms the existence of a fractional x, increasing with c, for various critical BEs. Two such cases are shown 
in figure 4, with their x~ results in close agreement with our analytical prediction, namely, eqs. (35, 41). 

VI. CONNECTION WITH PRBM MODEL 

In past, a random matrix ensemble, namely, power law random banded matrix (PRBM) ensemble was suggested 
as a possible model for the critical level statistics of Anderson Hamiltonian [6] . A PRBM ensemble is defined as the 
ensemble of random Hermitian matrices with matrix elements Hij as independently distributed Gaussian variables 

with zero mean i.e < >= and the variance < Hfj. s >= G^ 1 [l + (\i — j\/b) 2 ~\ , GV, = f3(2 — 5{j) and Gij = 1/2. 
It is critical at arbitrary values of the parameter b and is believed to show all the key features of the Anderson critical 
point, including multifractality of eigenfunctions and the fractional spectral compressibility. 

The success of PRBM ensemble, a one dimensional system, as a model for Anderson systems in arbitrary dimension 
is a little surprising. However, it can be explained on the basis of our formulation. The PRBM-AE connection is a 
special case of our study connecting any generalized Gaussian ensemble with BE. The PRBM ensemble being Gaussian 
in nature, its complexity parameter can be defined by using eq.(4) which can then be used to obtain its BE analog. 
The eq.(4) gives (with 7 = /3) 

Y - Y °= N{N U-(3) t} N - rMl + {bM21 

(42) 

which gives Y - Y a oc f(b)/N with f(b) w 26 85 ln(56) for b » 1 and f(b) w 26 175 for b « 1. As Y - Y w o(l/N), 
the R\ can then be given by eq.(14). Following eq.(15), A for a PRBM ensemble is 
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A prbm (b,E) = f(b)F 2 {E)C, 2 N-\ 



(43) 



The well-known size independence of level-statistics for PRBM case for all 6 values requires A pr t, m to be N indepen- 
dent which gives £ oc N 1 / 2 for all b ranges by our formulation. However note that, for PRBM model, < I 2 >oc N~ D2 
with D 2 as a function of b [18,19,35]. The use of ( oc< I 2 > _1 in eq.(43), therefore, gives a size-dependent A pr i, m . 
Keeping in view the well-known criticality of PRBM system for all b ranges, it seems that the relation ( oc< I 2 > _1 
is not valid for the PRBM case. 

Using the prediction ( — CqN 1 / 2 in cq.(43) and the relation A pr b m = A(,, the level-statistics of a PRBM can be 
mapped to a critical BE ensemble with 

c = (47r/3A prbm )- 1 c- £2 = (4^/(6)C 2 ^ 2 (^)e £2 )- 1 (44) 

The spectral statistics of PRBM therefore shows a crossover from from Poisson (as c — > oc for Ap r ^ m — > i.e b -> 0) 
to Wigncr-Dyson behavior (c — > for A pr i, m — > 00 or b — > 00). 

The spectral compressibility x for a PRBM ensemble at E — can now be obtained by substituting A = A pr b m = 
f(b)Co in the eqs. (34,35,41) which give 

X = 1 - 47T 2 C 2 /(&) for 6 « 1 

= (/ftr^/Wr 1 for6»l. (45) 

The above results are at least in the same form as obtained in [18,19,35]; the lack of explicit knowledge of £0 prevents 
us from making any further comparison. As \, in the eq.(45), changes from to 1 with decreasing 6, it violates the 
relation (27) in the range b « 1. A same violation was observed in previous PRBM studies [18] too. Thus our results 
obtained by using PRBM-BE analogy seem to be in accordance with earlier studies on PRBM model. 

In brief, the PRBM ensemble, with b as a parameter, can be mapped to critical BE with parameter c (see eq.(44)). 
As a consequence, the studies suggesting the analogy of spectral statistics for the PRBM and AE ensembles are well in 
agreement with our study claiming the AE-BE analogy. By using the connection of PRBM with MNS model [29,34], 
the PRBM-BE- AE analogy can further be extended to MNS-BE-AE analogy. 



VII. CONCLUSION 



In the end, we re-emphasize our main result: 

Under independent electron approximation, the level-statistics for the disordered systems undergoing localization — ► 
derealization transition of wavefunctions can be described by the Brownian ensembles (with uncorrelated elements) 
undergoing a similar transition. 

The analogy helps us in making following deductions: 

(i) The transition in the statistics is governed by a single scaling parameter A = ( j;) 2d = /(£)■ The second 
equality follows from the dependence of wavefunction statistics e.g inverse participation ratio I 2 and therefore £ on 
the complexity parameter | Y — Yq | . 

(ii) The level-statistics is governed by the competition between complexity parameter and local mean level spacing. 
The critical point of level statistics occurs when the complexity parameter Y — Yq and A; have same size dependence. 
In particular, if \Y— Yq\ <~ o(N a ) and A; <~ o(N^) for a disordered system, its critical point will occur when a — 2(3 = 0. 
However if the local mean level spacing in the system changes at a faster rate with size as compared to \/Y — Yq (i.e 
> a/2), the system will never reach its critical state and will always remain in the localized regime. 

(iii) The critical BE analog of a critical AE is unique. Further, it is different for critical AE with and without 
time-reversal symmetry. Similar to AEs, the level-statistics of BEs shows a scaling behavior as well as a critical point 
with fractional level-compressibility. However, unlike AEs, the x turns out to be symmetry dependent for BEs, their 
parameter A being symmetry independent. 

(iv) The AE-BE analogy confirms the symmetry-independence of the compressibility of levels and the multifractality 
of the wavefunctions at the critical point of Anderson transition. The analogy also indicates the non- validity of relation 
D 2 = d(l — 2\) in the strong multifractality limit, and, the approximate nature of the relation k = (2%) _1 . 
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(v) The AE-BE analogy helps us in formulating, for the first time, the exact 2-point level density correlation at the 
critical point of a disordered system. The formulation is applicable for a wide range of system parametric conditions. 

It should be noted that both MIT as well as WDT occur due to derealization of the wavefunctions. In fact, 
our analytical study suggests that the level-statistics of almost all complex systems undergoing a localization — > 
derealization transition follows the same route although with different transition rates; the state of level-statistics 
of two systems with different complexity may correspond to two different points on this route. In principle, our 
analytical work is applicable to the Gaussian models of complex systems only, however the intuition based on earlier 
studies suggest the validity of the results for the systems with other origins of randomness too [3] . For example, the 
investigation of a number of dynamical systems seems to support this intuition. It has been shown that the spectral 
statistics of pseudo-integrable billiards is remarkably similar to the critical statistics of AE [37]. The presence of 
a statistics intermediate between the Poisson and GOE has already been shown for the Kicked rotor in the non- 
integrable regime of the kicking parameter [38] . A correspondence of the integrable systems to the insulators and of 
the chaotic systems to the metals is already known to exist. The Integrability — > chaos transition in the dynamical 
systems therefore seems to follow a route in the level-statistics similar to that of MIT; note such a transition in 
classical systems corresponds to a derealization of the wavefunctions in their quantum analog. Thus the analogy of 
the statistical level fluctuations between AE and BE may possibly be extended to dynamical systems and BE too; if 
the latter is found to be correct, the analogy would be useful to obtain the correlations for the non-integrable regime. 

The evidence of such an analogy would suggest the existence of several features, unknown so far, for the level 
statistics of dynamical systems. For example, the analogy can be used to intuitively claim and search for the existence 
of a critical point, the dimensional dependence of level-statistics and the multifractality of eigenfunctions during 
the transition from Integrable to chaotic dynamics. It should be noted that a generic one-dimensional dynamical 
system always shows a Poisson level-statistics (in analogy with one-dimensional AE). However the dynamics in a 
3-dimensional system shows a feature namely "Arnold Diffusion", absent in lower dimensions. The intuition based on 
the above analogy suggests the possible existence of a critical level statistics at the parametric values at which Arnold 
Diffusion takes place. A further exploration of such an analogy is therefore highly desirable. 
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FIGURE CAPTIONS 

Fig. 1. The study of critical BEs: 

(a) The scaling behavior of the integrated nearest-neighbor spacing distribution a during WDT. Note a values for 
BEs with different parameters /i and sizes N converge on the same two curves, thus indicating a dependence on a 
specific combination of fj, and N, namely, z = n/N 2 . Further, at z — c, a remains unchanged for different N values, 
thus indicating a critical point of BEs. 

(b) The comparison of tail of P(S) distribution for two of the critical BEs with function c~ KS for (3 = 2 case. For 
intermediate c-ranges, the P(s) is well-fitted by the function c~ KS . The fitting however, seems to be poor for smaller 
c values which is as expected, due to statistics approaching GUE limit (which corresponds to P(s) <~ c _7rs ' 4 ). The 
fitted K-values are as follows: — 0.8 for c = 0.3, (ii) n = 1.7 for c = 0.03. The above K-values seem to deviate 
significantly from the relation k = (ir/A-yc) (obtained by using k = (2x) _1 with \ given by eq.(35) and 7 = 2). We 
have seen a similar deviation for the BEs with (3 = 1 too. This suggests the non- validity of relation k = (2x) _1 in 
general although it seems to work for some c values (see for example figure 3). 

Fig. 2. The behavior of level-density F(E) = N' 1 .^: 

(a) for the two cases for two system sizes L = 10 and L = 13. The numerical fitted function has the form 
F = hc- {E2 /f^ with /1 = 0.16, f 2 = 5 for AE t and /1 = 0.016, f 2 = 400 for AE„ t . 

(b) for the critical BE analogs of the higher-order correlations of the two cases considered in 2. (a). Here F for 
all the critical BE cases is well-fitted by the function F(E) = tt^ 1 / 2 c^ e . Note the lack of analogy between the 
mean-level densities for the cases given in 2(a) and 2(b) while their higher order correlations (shown in Figures 3,4) 
are approximately same. 

Fig. 3. The comparison of distribution P(S) of the nearest-neighbor spacings S for the AE (d=3) and BE cases on 
a log-log scale. To confirm the critical state of the AE, the distribution is shown for two system-sizes L for each AE 
case. The insets show the same functions on a lin-log scale and also compares the behavior with e~ KS : 
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(a) AE t (with hard wall boundary conditions, random hopping and time-reversal symmetry) and its critical BE 
analog (c = 0.1). The dashed line in the inset is fitted function / = 4c~ 1 - 7S which gives n 1.7. (b) AE nt (with 
periodic boundary conditions, non-random hopping and no time-reversal symmetry) along with its critical BE analog 
c = 0.2. The dashed line in the insert is fitted function / = 2e _1 which gives k w 1.5. 

Fig. 4. The comparison of the E 2 (r)/r-behavior for the AE and BE cases: (a) AE t and the corresponding critical BE 
analog (c = 0.1). (b) AE nt and corresponding critical BE (c = 0.2) 

Note here the critical BE analog for each AE case is same as for the P(5')-study. As can be seen, E 2 /r for large 

r seems to approach the limit suggested by the relation \ = ^^oo^ 2 ^ = ^ C /^ 7T ^ is, x — 0-25. Note this 
is the expected x f° r AEs on the basis of eqs. (25,27) too. Besides showing AE-BE analogy, the figure also confirms 
that (i) symmetry independence of \ f° r AEs (ii) fractional x~ result for a critical BE. 
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